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Abstract 

We investigate the dynamics of a macroscopic system which consists of an anhar- 
monic subsystem embedded in an arbitrary harmonic lattice, including quenched 
disorder. The coupling between both parts is bilinear. Elimination of the harmonic 
degrees of freedom leads to a nonlinear Langevin equation with memory kernels 
r(t) and noise term C,{t) for the anharmonic coordinates q(t) = {qa{t)). For zero 
temperature, i.e for C,{t) = 0, we prove that the support of the Fourier transform of 
r(t) and of the time averaged velocity-velocity correlations functions K[t) of the 
anharmonic system can not overlap. As a consequence, the asymptotic solutions 
can be constant, periodic, quasiperiodic or almost periodic, and possibly weakly 
chaotic. For a sinusoidal trajectory q{t) with frequency we find that the energy 
Et transferred to the harmonic system up to time T is proportional to T". If 17 
equals one of the phonon frequencies ojy, \i \s a = 2. We prove that there is a zero 
measure set C such that for 17 in its full measure complement TZ\C, it is a = 0, 
i.e. there is no energy dissipation. Under certain conditions C contains a subset 
C such that for Q £ C the dissipation rate is nonzero and may be subdissipative 
(0 < a < 1) or superdissipative (1 < a < 2), compared to ordinary dissipation 
(a = 1). Consequently, the harmonic bath does act as an anomalous thermostat, in 
variance with the common belief that elimination of a macroscopically large number 
of degrees of freedom always generates dissipation, forcing convergence to equilib- 
rium. Intraband discrete breathers are such solutions which do not relax. We prove 
for arbitrary anharmonicity and small but finite coupling that intraband discrete 
breathers with frequency 0, exist for all Q in a Cantor set C{k) of finite Lebesgue 
measure. This is achieved by estimating the contribution of small denominators ap- 
pearing for G{t;i}), related to T(t). For Q G C{k) the small denominators do not 
lead to divergencies such that G{t; 0) is a smooth and bounded function in t. 
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1 Introduction 



The analytical treatment of the dynamics of a macroscopic system of N parti- 
cles is not possible, in general. But there are situations which allow to simplify 
the dynamical description. We will restrict ourselves to classical systems. One 
of the prominent examples is Brownian motion. A big particle with mass M 
interacts with a solvent. If the solvent particles have a mass m much smaller 
than M a time scale separation exists. Separating the fast motion of the sol- 
vent particles from that of the big one, called Brownian particle, a Langevin 
equation results for the motion of the Brownian particle. The influence of 
the solvent particles occurs through a friction and a fluctuating force. These 
two are related to each other by the so-called fluctuation- dissipation theorem 
of second kind. For a general system with time scale separation, the Mori- 
Zwanzig projection formalism [1,2] allows to derive a Langevin equation for 
the slow variables. If there is a single slow variable A, then one gets [1,2] 




where f2 is a microscopic frequency, M{t) a memory kernel and f{t) the fluc- 
tuating force. Elimination of a macroscopic number of fast degrees of freedom 
always leads to memory effects. If M{t) relaxes much faster then the average 
of A{t), one can perform a Markov approximation, i.e. M{t) ^ 2^S{t — t'). The 
analytical calculation of M(t) is not feasible, in general. However, there exists 
a class of models for which the exact elimination of a macroscopic number of 
degrees of freedom is possible, even if there is no time scale separation. These 
models we call Kac-Zwanzig models [3,4]. They consist of a subsystem coupled 
to a macroscopic bath of harmonic oscillators. The subsystem can be micro- 
scopic, i.e. the number of its degrees of freedom is of order one or it can also be 
macroscopic with, e.g. 10^^ degrees of freedom. Since the equations of motion 
are linear in the oscillator coordinates, the harmonic degrees of freedom can 
be microscopic, i.e. the number of its degrees of freedom is of order one or it 
can also be macroscopic with, e.g. 10^^ degrees of freedom. As a result, one 
obtains for the coordinates of the subsystem a stochastic differential equation 
of type of eq.(l). The memory kernel M (which is a matrix in general) can 
be expressed by an integral over the frequencies of the harmonic bath. Its 
spectral properties strongly influence the time dependence of M and therefore 
that of the coordinates of the subsystem, too. In the Langevin case where the 
spectrum is absolutely continuous with infinite support, the subsystem will 
relax to equilibrium for all of its initial conditions. In this case, the harmonic 
oscillators act as a normal thermostat. However, if, e.g. that spectrum contains 
a dense discrete component as it occurs for harmonic systems with quenched 
disorder, the situation may be different. Then, there exist initial conditions. 
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which will not relax to equilibrium, but may converge to periodic, quasiperi- 
odic or possibly weakly chaotic solutions. This behavior is closely related to 
the existence of spontaneous energy localization. 

Spontaneous energy localization in large nonlinear discrete systems is now 
known as an ubiquitous phenomena related to the existence of Discrete Breathers 
(DBs). When a sufficiently large amount of energy is injected locally in a sys- 
tem which may sustain DBs and which is initially at low temperature, it is 
found that although some part of this energy spreads over the system, a sub- 
stantial part of this energy may remain localized as a DB [5] over very long 
time. The discovery of this kind of localized modes was so surprising because 
they existed in systems with discrete translational invariance in space. DBs 
also called Intrinsic Localized Modes (ILM) are spatially localized time peri- 
odic solutions of discrete nonlinear Hamiltonian systems [6,7,8,9]. Since many 
of them are linearly stable and thus can trap energy over long times, they 
should play an essential role in energy relaxation in complex systems. 

The existence of DBs has been rigorously proven by several methods and for 
a variety of models with optical or acoustic phonons, with or without gaps, 
in one or any dimension and without or with randomness [10,11,12,13,14]. All 
these existence proofs require that the DB frequency and its harmonics do not 
belong to the phonon spectrum (obtained by linearization of the dynamical 
equations in the vicinity of the ground state). Otherwise, the DB would be 
expected to radiate energy in the phonon band and to decay. In addition, 
the proofs [10,11,12,13,14] only concern the existence of out-of-band DBs, i.e. 
their frequencies and those of the higher harmonics are in the band gaps of the 
phonon spectrum. For certain models the existence of DBs was even proven 
before by Albancsc and Frohlich [15,16] and by Albancse, Frohlich and Spencer 
[17] under the requirement of randomness. In that case the eigenstates of the 
linear eigenvalue problem are already localized (in dimensions larger than two 
the randomness has to be large enough) and therefore it is more intuitive that 
localized time periodic solutions could persist in the presence of nonlinearities, 
although to prove this is highly nontrivial 

When the phonon spectrum is discrete, for example in a strongly disordered 
system with Anderson localization, energy radiation of linear phonons by a 
localized time periodic solution cannot occur. Then, the existence of Intraband 
Discrete Breathers (IDBs) with their frequency or some of their harmonics 
inside the phonon spectrum, might be a priori possible. Indeed, numerical 
evidence [18,19] suggested that such solutions do exist. More precisely, it was 
conjectured that for each localized Anderson mode, there is an associated 
family of IDBs with variable frequencies (and amplitude), which converges to 
this Anderson mode at the small amplitude limit. It was also conjectured that 
the IDBs with respect to their frequency dependence do not form a continuous 
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family of sohitions as for ordinary (out-of-band) DBs, but are defined only on 
a fat Cantor set of frequencies (i.e. with nonvanishing Lcbcsgue measure). A 
forbidden frequency gap is associated with each of the other frequencies of 
the linear Anderson modes. These frequencies are dense in the linear phonon 
spectrum but the width of these gaps drops to zero very fast as a function 
of the spatial distance of the corresponding Anderson mode of the center to 
the IDB. Then, the complementary Cantor set of frequencies where the IDB 
exists is non void with a strictly positive Lebcsguc measure. Moreover, the 
measure density of this Cantor set is expected to go to full measure when the 
amphtude of the IDB goes to zero that is when the IDB tends to be a linear 
mode. 

The IDB should not be confused with so-called embedded solitons in discrete 
lattice systems (see [20,21] and references therein). These are nonlinear exci- 
tations with a discrete frequency lying within the linear spectrum, whereas 
IDB take frequencies from a set of full measure. Furthermore, in contrast to 
IDB the existence of embedded solitons requires exponentially localized linear 
modes. 

The first proof of existence of such IDBs for a certain class of random mod- 
els was done by Albanese and Frohlich [15,16] and by Albanese, Frohlich and 
Spencer [17]. The models studied by these authors are the random, discrete 
Schrodinger equation and the wave equation with a cubic nonlinearity. A per- 
turbational approach was used to prove that all periodic solutions of the linear 
equation can be continued for small enough A, the strength of nonlinearity. 
For arbitrary A it was proven that there exist at least certain periodic solu- 
tions. Transforming a random, nonlinear dynamical system without acoustic 
phonons into Anderson space, i.e. to the space of the localized eigenstates of 
the linear problem and restricting to nearest neighbor interactions in Ander- 
son space, Frohlich, Spencer and Wayne [22] have even proven (under certain 
conditions) the existence of localized almost periodic solutions, i.e. of KAM 
tori. 

Besides the existence of IDBs, it was also conjectured in [18,19] the existence of 
families of multibreather states which exhibit several energy peaks associated 
with different Anderson modes. We may explain intuitively this conjecture by 
the fact that the frequency of a nonlinear mode depends on its amplitude. 
It is thus possible to tune the amplitude of an arbitrary number of differ- 
ent Anderson modes in order they get the same frequency. Despite nonlinear 
modes cannot be simply superposed as in the linear case, a time periodic 
multibreather solution can be constructed numerically. When, the number of 
IDBs involved in the multibreather state is infinite, these multibreather states 
are spatially extended and able to transport energy by phase torsion [23]. As 
a result, nonlinearity restores the ability of the system to transport energy. 
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The aim of this paper is to provide an exact proof of these conjectures for 
a special class of random models where the underlying physical arguments 
appear explicitly. These models, consisting of a microscopic or macroscopic 
anharmonic subsystem coupled to a bath of harmonic oscillators, are designed 
on one hand for obtaining simpler mathematical proofs and on the other hand 
to contain nevertheless the essential features, e.g. for the trajectories q{t) of 
the anharmonic system, which can be found in more general complex and re- 
alistic models. Although these results do not provide a rigorous proof (in the 
mathematical sense) for general and fully anharmonic systems, they neverthe- 
less make clearer the fundamental reasons for their existence and thus support 
the earlier conjecture that they might exist in general. 



Particulary, we will derive for rather general anharmonic and harmonic subsys- 
tems(arbitrary dimension, with or without disorder) an exact formula which 

relates the energy dissipation rate Et of the anharmonic system averaged over 
a time interval of length T to the Fourier transform r(a;) of the dissipation 
kernels T{t) and to K{uj), the Fourier transform of the velocity-velocity corre- 
lation functions of the trajectories q{t). For finite initial energy and bounded 

Hamiltonian from below it is Et — > for T — > oo. This will yield the nonover- 

lapping criterion for the support of T{uj) and K{uj) which will have conse- 
quences for the possible asymptotic behavior of q{t) for t ^ oo. Calculation 

of the energy dissipation rate Et for a sinusoidal force with frequency Q 
will demonstrate the existence of extreme sensitivity of the dissipated energy 
by that force on its frequency ft. Compared to ordinary dissipation where 
Et ~ T, superdissipative and subdissipative behavior may occur for which 
Et increases, respectively, faster and slower than T. We will also show that 
there exists a set of full measure such that Et is bounded for all Q in that set, 
i.e. no energy dissipation exists at all. 



The possibility of no energy dissipation allows the existence of periodic solu- 
tions q{t + 27r/il) = q{t), also within the phonon band. We will show that 
such IDB exist with frequencies on a fat Cantor set. In that case relaxation to 
equilibrium is prohibited such that the harmonic bath acts as an anomalous 
thermostat. This behavior is related to the pioneering work of Fermi, Pasta 
and Ulam, which has given evidence for nonergodic dynamics. New results 
confirm that ergodicity in phase space of large and complex dynamical sys- 
tems is not as smooth as it was believed through Langevin bath theory. When 
starting from initial states which are not at thermal equilibrium, there are 
many islands in phase space which may trap the state of the system for long 
time and could result in unusual and slow relaxation processes toward the 
thermal equilibrium. 
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2 Modeling the Intraband Discrete Breather Problem 

In order to keep our presentation selfcontained we review and discuss prop- 
erties of linear spectra which have been proven earlier or are well-admitted 
in order to define a model and an associated set of hypothesis which could 
match with the different physical situations. Our purpose is to describe the 
main characteristics which could be expected for large families of nonlinear 
random models with linear spectrum involving both discrete and continu- 
ous components. Continuous spectra occur only in the thermodynamic limit, 
i.e. for an infinite system. 

In order to fix the ideas, let us consider as an example a random model on 
a d-dimensional hypercubic lattice, for example a Klein-Gordon type model 
with Hamiltonian 

^ = E (f + + E ' (2) 

where the momentum pi is the variable conjugate to the displacement Ui at 
site i. The local potential Vi{ui) is random. For example it could have the form 

V^{u,)^^ul + ^uU^< + ... (3) 

where uJi is chosen according to some probability law P{uj). denotes the 
bonds which connect neighboring sites (i, j) (counted once). Accordingly, the 
sum in eq. (2) is a double sum over all i and j which are nearest neighbors. C is 
the constant of the harmonic coupling between nearest neighbors, which is as- 
sumed to be identical for all nearest neighbor bonds. The dynamical equations 
are 

Ui + V:{ui) + J2C{u,-Ui) = Q (4) 

j:i 

j : i denotes the sites j which are nearest neighbors to i, i.e. the sum in eq.(4) 
is a single sum over j which are nearest neighbors of i. DBs are spatially 
localized time periodic solutions of eq.(4). Thus, their amplitude must decay 
to zero at infinity. Consequently, we may approximate this dynamical system 
by another dynamical one which obeys the same eq.(4) for the few sites i 
where the amplitudes of the DB is not small and the linearized equations 

Ui + Vf{Q)ui + J2 C{uj -Ui) = Q (5) 

j:i 

for all the other sites i. Thus, the system becomes equivalent to a finite an- 
harmonic system coupled linearly by local interactions to a harmonic system 
(phonon bath). Note, we start with a finite harmonic system and finally take 
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the thermodynamic hmit. For a finite harmonic system the hnear spectrum is 
always discrete. 

This kind of approximation strongly simphfies the investigation of DBs. As- 
suming the existence of time periodic solutions (DBs) of the isolated anhar- 
monic system (i.e. uncoupled with the harmonic phonon bath) we will discuss 
under which conditions these solutions could persist when the coupling to the 
phonon bath is non vanishing, although not too large. This harmonic phonon 
bath considered uncoupled to the anharmonic system, is characterized by its 
frequency spectrum. When the frequency of the DB and its harmonics do not 
belong to the frequency spectrum of the phonon bath, the proof that the DB 
solutions can be continued to nonzero coupling can be obtained as in [10]. 
This will be performed in the 5. section. When the frequency of the DB, one 
of its harmonics or several belong to the frequency spectrum of the phonon 
bath, it is essential to consider whether this phonon spectrum is discrete or 
continuous. 

Figs. 1 and 2, respectively, show possible schemes representing situations where 
the hnear phonons are either spatially extended or localized at the DB fre- 
quency. 




Fig. 1. The eigenmodes of the uncoupled phonon bath are extended waves at the DB 
frequency. Then, at nonvanishing coupling periodic oscillations of the anharmonic 

system should generate waves which propagate and carry energy at infinity. As a 
result, the DB solution cannot persist as a finite energy (spatially localized) time 
periodic solution. 

The model we propose is simply to consider the system as harmonic in the 
regions where the amplitude of the DB is expected to be small and to keep it 
anharmonic in the region where the solution has a large amplitude. Of course, 
our model shall depend on the location of the IDB or multiDBs we consider 
since it determines the selection of the anharmonic subsystem. 

This kind of simplification also allows to investigate the relaxation of an arbi- 
trary localized initial excitation. If the finite initial energy is spread completely 
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Fig. 2. The eigenmodes of the uncoupled phonon bath are spatially localized at 
the DB frequency. No energy propagation is possible. Then, the DB solution of the 
anharmonic system may persist at finite coupling. 

over the system the displacements become arbitrary small and the lineariza- 
tion of eq.(4) is allowed for large enough times. However, if the spreading is 
incomplete, we are allowed to linearize eq.(4) outside the localized region. 

When {ui} is uniformly small, the linearized equations which yield the linear 
phonon spectrum are 

€i + u;^ei + J2C{ej-ei)^0 (6) 

j:i 

which has time periodic solutions of the form e\°'\t) = e-°^e*'^"*. The eigenfre- 
quencies and eigenmodes 6^°^ = {^1"^} obey the eigenequation 

(Aoe("))^ = u^eT^ + ^ ^(ef - e^^) = c^ef) (7) 

j:i 

where Aq is the matrix of the second order variation of the potential energy 
from its ground state. The eigenfrequencies uja should not be confused with 
the frequencies cui of the small amplitude oscillations in the local potential 

2.1 Linear Spectrum Properties 

It has been proven [24] that the linear spectrum of eq.(7) is discrete and its 
eigenvectors (Anderson modes) are square summable and exponentially local- 
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ized with probability one, when either the spatial dimension is less than three 
or the coupling is not too large. At larger coupling and when the dimension 
of the lattice is at least three, it is believed that there exist mobility edges 
in the spectrum (although no rigorous proof is available up to now). There 
are two mobility edges which determine a central interval where the spectrum 
is continuous. The corresponding eigenmodes are extended and non square 
summable. In two intervals limited by the edge of the spectrum, the spec- 
trum is discrete and the corresponding eigenmodes are square summable and 
exponentially localized. 

More generally, we may consider nonlinear models on arbitrary discrete lattices 
where the small amplitude solution leads to the eigenvalue equation 

Ax^ = EaX^ (8) 

for an arbitrary positive bounded self-adjoint linear operator A operating in 
the Hilbert space x = {xn} where n represent the sites of an infinite discrete 
lattice C in 1,2,3... dimensions (or any discrete network), with the standard 
Hermitian product (x|y) = J^n^nUn and ||x|| = ((x|x))^/^. Since we consider 
physical systems which have to be linearly stable, we must assume the spec- 
trum to be positive or zero and then, the eigenvalues Ea = uj"^ are the square 
of the eigenmode frequencies. We may also consider systems with acoustic 
phonons where the spectrum is gapless and contains the translation mode 
Xn = a associated with zero eigenfrequency which is always extended. 

Standard mathematical theories decompose the spectrum S of such a self- 
adjoint bounded linear operator A in three parts: the absolutely continuous 
spectrum, the discrete spectrum and the singular continuous spectrum [25]. 
The characteristics of the spectrum manifests itself in finite size systems in 
the limit of large size N as explained now. The eigenvalue equation restricted 
to a finite subset Af of N indices of C corresponds to the diagonalization of a 
finite self-adjoint matrix 

A(^)x(^) = E(^)x(^) (9) 

which yields orthogonal and normalized eigensolutions x^^ with compact 
support A/" corresponding to real eigenvalues. Taking A^ —>■ +oo such that 
the finite support tends to cover the infinite lattice C, the nonvan- 

ishing accumulation points for the weak topology of the eigensolutions 
of finite systems x^^ are eigensolutions of eq.(8) for an eigenvalue Ea, the 
corresponding limit of E^\ The nonvanishing limit solutions are necessarily 
square summable with a nonzero norm smaller or equal to ||x|;^)|| = 1 and 
thus can be normalized. 

^ A vector is said to converge to a limit vector for the weak topology when each of 
its components converge to a limit. 
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This set of square summable eigensolutions with eigenvahies determines 
the discrete spectrum of A. Since C is countable, the discrete set of normahzed 
eigenvectors is also countable (or finite). The inverse rate of the exponential 
decay of the eigenmodes determines the so called localization length which 
is a smooth function of the eigenfrequency and diverges when approaching 
mobility edges. In systems with acoustic phonons in 3 dimensions and more, 
it is believed that the spectrum becomes continuous in a nonvanishing positive 
interval close to zero and then there is only one mobility edge. 

Extended states (e.g. plane waves) are not square summable and their norm 
diverges proportionally to the system size N. These eigenstates can be nev- 
ertheless obtained as nonvanishing accumulations points of \/iVx^^ This set 
of nonsquare summable solutions and their corresponding eigenvalues of the 
eigenvalue equation (8) are usually called pseudo eigensolutions and pseudo 
eigenvalues. They corresponds to the absolutely continuous spectrum of A. 
They form uncountable continua. 

The limits of the other finite size eigenstates which behave spatially differ- 
ently for example which develop a fractal structure at large scale or which 
exhibit algebraic decay at infinity, determine the remaining part of the spec- 
trum which is the singular continuous spectrum of A. Although there are 
examples of operators with purely singular continuous spectrum (for exam- 
ple with quasiperiodic potentials), their contribution to the spectrum of most 
operators is expected to be marginal, generally, and then it reduces to few 
isolated specific points of the spectrum (corresponding to the mobility edges) . 

The spectrum of an operator A (for example random) is a closed subset S 
of the real axis. It might be purely absolutely continuous (for example for 
operators with spatial periodicity) or purely discrete (for example for operators 
in Id with randomness or higher dimension with strong disorder) or it might 
decompose into two disconnected subsets Sd and Sac where the spectrum is 
discrete and absohitely continuous, respectively, separated by mobility edges. 
This situation is supposed to occur for random operators in 3 dimensions and 
more at weak enough disorder [19]. 



2.2 The anticontinuous Limit 

Considering now again, the nonlinear system (4) at the anticontinuous limit 
C = 0, it simply becomes a collection of uncoupled nonlinear oscillators with 
Hamiltonian Hi = ^ + Vi{ui) and thus is trivially integrable. It is then conve- 
nient to represent each anharmonic oscillator with the standard action-angle 
variables Ii,6i. Then, the local Hamiltonian Hi[Ij) is only a function of its 
action 7^, and its frequency is Hl{Ii) = Ui{Ii). Note, that uii{Ii) should not be 
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confused with the frequencies Ui in eq.(3). Time periodic solutions at frequency 
VL are obtained by choosing each oscillator i either at rest or oscillating at fre- 
quency Vt = uJi{Ii) if there exists /j such that this equality may be fulfilled. 
The solution corresponding to a Discrete Breather at site i and at frequency 
Q is obtained by choosing all the oscillators j i a,t rest, except oscillator i 
which is oscillating at frequency 

At C = 0, the linear spectrum of eigenfrequencies defined by eq.(7) is the 
closure Sl = Ui{'-^i} of the set of linear frequencies. Considering a DB solution 
at the anticontinuous limit at site i and at a frequency which is nonresonant 
with the linear spectrum (that is pVt ^ Sl for any integer p), and fulfills the 
condition that duji{li) / dli = H'l{Ii) 7^ at uii{Ii) = Q, it was proven by 
the implicit function theorem [10], that this solution can be continued as a 
spatially localized solution for C 7^ belonging to a nonvanishing interval. 

When Q G Sl, the proof for the DB existence given in [10] which requires 
the invertibility of Aq — Q^l breaks down. However, numerical investigations 
suggested the existence of stable IDBs with frequency fl E Sl although they 
are not continuous functions of the model parameters. Otherwise, this numer- 
ical work also suggested the existence of intraband multisite breathers called 
multibreathers. A direct proof of existence of IDB in fully anharmonic mod- 
els is likely possible (as an extension of the present work) but it should be 
more complex for technical reasons than in the simplified model we propose. 
We think, however, that this model captures the essential physical features of 
general nonlinear random models. 



2.3 Defining a model for IDBs 

We denote by q = the set of variables which describe the coupled an- 
harmonic oscillators and by p their conjugate variables. This subsystem is 
described by an anharmonic Hamiltonian 



Without restricting generality the masses are chosen to be equal to one. The 
single anharmonic oscillator where q is a scalar, is the simplest but a special 
case. This anharmonic part is linearly coupled to an infinite harmonic system 
described by the set of variables u = {ui} and conjugate momenta p = {pi} 
where the indices i are nodes of an arbitrary network. The harmonic Hamil- 
tonian is HhaTTaiiuijPi}) = \p'^ + ^u^Mu where M is a bounded self-adjoint 
operator and t denotes transposition. Again, the masses can be taken equal 




(10) 
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to unity without lack of generality. The total Hamiltonian takes the form 

H = Hoip, q) + Hharm({Mi,]?^}) " Aq*Cn (11) 

where the matrix of coupling constants is C = {Ca,i} and A is introducted for 
convenience. These coupling constants are short ranged that is J2i |C'a,jp < 
+00. 

It is more convenient to describe the harmonic system as a collection of nonin- 
teracting harmonic oscillators u which corresponds to its normal modes. The 
normal mode coordinates of the harmonic oscillators of the phonon bath are 
now the coefficients u^, of the expansion of {ui} with respect to the normalized 
eigenmodes u^''^ with frequency lOi, which fulfills Mu^"^ — ujIu^^\ They are 
now coupled to the anharmonic system by the set of coefficients Cj, = u^'^^^C 
or explicitly Cj.,a = T,u\''^Ca,i- 

i 

As we assumed above, the set of coupling constants Ci = (C«,j) between q 
and Ui are short ranged. In most physical situations, they decay exponentially, 
that is there exists two positive constants K and 7 such that 

|Ci| < i^e-^l'l . (12) 

Thus, we may assume the stronger Li condition |Ci| < +00 without loosing 
much generality. 

Then, we obtain a total Hamiltonian of the form 

7^ = 7^o(p, q) - A E Ctq u^ + T, (liil + ^c^X) ■ (13) 



The reorganization energy Hr is the change in the groundstate energy of the 
whole system under the constraint of a static displacement q of the anharmonic 
oscillators. It is obtained by minimizing the last two terms of the energy (13) 
with respect to u^, for fixed q. This readily yields Hr{(\) — —^q*Kijq where 

In well-behaved physical models, this sum is always finite since any finite 
static displacement q should involve definite reorganization energy (and a finite 
force) . This condition requires the assumption that the positive definite series 
((14)) is always absolutely convergent. For example, the case of models with 
acoustic modes where one of the uOy may vanish because there is always a mode 
at zero frequency (which corresponds to the free motion of the center of mass of 
the whole system) . This feature may cause a divergency for coefficients of series 
(14). But in the limit of low frequency, these (acoustic) phonon modes decouple 
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from the local degree of freedom describing the anharmonic part of the system 
(lim^^_^o C,^ = 0) and thus this reorganization energy (14) remains finite since 
Cu is linear or even of higher order in uju for ujn 0. The reorganization 
energy Tiji renormalizes the potential energy of the anharmonic subsystem, 
i.e. we can introduce an effective energy 7Ye//(p, q) = 'Wo(p, Thus, 
the total Hamiltonian (13) may be written as 

7^ = K//(p,q) + E(^^^ + ^^^K-^f) , (15) 

where Uy{c\) = AC^q/a;^ is the harmonic displacement minimizing the last two 
terms of (13), for fixed q. Since we did not specify properties of the harmonic 
subsystem, which is considered as very general, it is now necessary to introduce 
some assumptions at this level for ensuring our model to be physically realistic. 
The density of the positive eigenfrequencies uj„ of the harmonic bath (called 
density of states) is the limit of the positive measure defined as 

In physical situations, it is generally an absolutely continuous measure, dense 
on a finite interval (that is dniuS) — g{u!)du! where g{u!) is a measurable pos- 
itive function). But this assumption is not necessary and moreover does not 
characterize whether the phonon spectrum is discrete or continuous. 

For example, in spatially periodic harmonic systems the whole phonon spec- 
trum is absolutely continuous with extended eigenstates which are planewaves. 
In random Anderson models in three dimensions and at weak enough disor- 
der, the phonon spectrum S may split into two disconnected sets S = Sac U 
which consists of a central interval Sac where the spectrum is supposed to be 
absolutely continuous (corresponding to extended eigenstates) and the union 
Sd of two side intervals where the spectrum is supposed to be discrete (and 
corresponding to square summable eigenstates). These intervals are separated 
by mobility edges (isolated points) . At large disorder in 3d or at lower dimen- 
sion for any disorder, the whole spectrum is discrete and consists of a single 
interval Sa- 

Actually, in all these examples, the density of states is an absolutely contin- 
uous measure which does not determine a priori whether or not the phonon 
spectrum is continuous or discrete. Note also that in the case of random po- 
tentials, the density of states depends only on the statistical properties of 
the potential and not on the realizations. The characteristics of the spectrum 
(absolutely continuous, singular continuous or discrete) appears in our model 
through the properties of the distribution of coupling constants C^. They de- 
termine the properties of the kernel which appears after elimination of the 
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harmonic variables and which shall be discussed later. 



3 Extended Langevin Theory 

Models for an anharmonic system coupled to a harmonic phonon bath are 
ubiquitous in physics for describing relaxation processes, for example chemi- 
cal reactions within Kramers theory [4,26]. However, in those applications the 
spectrum of the harmonic bath was always assumed to be absolutely continu- 
ous. In that case, standard calculations yield after elimination of the phonon 
bath variables, extra damping terms for the anharmonic system (plus Langevin 
random forces, if the temperature is not zero). When there is a discrete spec- 
trum, the results are totally different, as we shall see now. 

The well-known advantage of a coupling with a harmonic phonon bath is that 
it is straightforward to eliminate the variables of the phonon bath in order 
to get effective equations for the anharmonic variables [4] (we already took 
advantage of such a situation for proving the existence of (interband) Discrete 
Breathers in systems with acoustic phonons [12]). 



3. 1 Derivation of a Langevin Equation 

Making use of eqs. (10) and (13), the dynamical equations are 



where V'{q} stands for the gradient of V with respect to q. The second equa- 
tion is linear and readily yields 




(16) 



V 



-\- ujIu^ - ACf,q = 



(17) 




(18) 



where 



u 



X 



-Ctq(O) + e^(0) cos uja + 



e.(0) 



sin oUi^t 



V 



v 
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is an arbitrary solution of the homogeneous part of eq.(17)[l]- If we assume the 
phonon bath in thermal equilibrium, e^^-* (t) is a random variable due to ther- 
mal fluctuations and we have at temperature T, {e^\0)el^\t)) = cosu^t. 
Substituting u„{t) from (18) into (16), we obtain 



q(t) + F'(q(t))-A2 / G(t-r)q(r)dr-A2^a(Ctq(0))cosa;,t-AC(t) = 



(19) 



where the coefficients Ga,i3{t) of matrix G(t) are defined by 

G^M = Y.^^^^^smuja (20) 



This series is absolutely convergent and yields a smooth function of t because 
of hypothesis (14) and < < oo: 

\n n I 

<^^E-J<+o« 

C(t) is the Langevin random force 

CW = E4°nt)c. (21) 



It is convenient to rewrite eq.(19) as 

t 

q(t) + K'//(qW) + A' / r(t - r)q(r) dr - K{t) = (22) 



where 



Veffiq) = V{q) - Y: ro.AO)qaqp (23) 

r«,/.W=E^^^cosa;.t (24) 

and dT{t)/dt = —G{t). We note r(0) = K/j which corresponds to the coeffi- 
cients of the reorganization energy (14). Eq.(22) has the form of the Langevin 
equation (1), however with the nonlinear term V^jj(qr(t)). 



^ Note that we choose the harmonic displacement €^{0) cos uj^t + ^^j^ sinujyt with 



respect to the equilibrium position -4C^q(0) for given initial coordinates q(0). 
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The damping kernel r{t) is related to the random force by the general 
Langevin relation 

(Ca(0)C/3(t)) = fc^T^ COSCU.t = kBTTg,f,{t) (25) 



where ( ) denotes the canonical average with respect to {£,^(0)} and {e,y{0)}. 

Definition (24) shows that the damping kernel T{t) = / r(ci;)e*'^*(ia; is the 
Fourier transform of r(a;) with matrix elements which are measures 

C C 

fa,/3(a;)da; = daa,p{uj) = tt^ "'"^"'^ {6(0; - u;„) + S{u; + u;t,))du; (26) 

1/ ^^11 



Actually, measure aa,p{u) may be written in a more compact way as the result 
of an integration in the complex plane 

a.A^) = [ da^A^') = < C.\ ^^^\^^ \C, > dz (27) 

on a closed contour C(a;) which encompasses the part of the spectrum of the 
positive linear phonon operator M on the real axis which is smaller than 
cu^. Here it is Cq, = {Ca,i)- Note that in the case of random potentials, this 
measure does depend on the realization of the disorder unlike the density 
of states. This positive measure daa,i3{(^) can be split into three components 
which are respectively absolutely continuous, discrete and singular continuous 
measures. eq.(27) proves that this decomposition precisely corresponds to the 
decomposition of the spectrum S of operator M which describes the phonon 
bath, into its absolutely continuous Sac, discrete Sd and singular continuous 
parts Ssc- Thus, the support of each of these measure components on the 
real cu axis is independent of a and (3. Consequently, the properties of the 
Fourier transform measure of the memory kernel T(ou) are directly related to 
the spectral properties of the phonon bath. 



3.2 Relation between Phonon Spectrum and Kernel Properties 



We now show how the properties of the kernel T{t) and of its Fourier trans- 
form r(a;) are related to the properties of the phonon bath. Actually, the 
Fourier transform eq.(26) of the kernel T{t) is not a priori a smooth function 
but is a positive measure T{uj)duj which very generally can be decomposed 
into three parts: the absolutely continuous part Tac{'-^) where Tad^^) is a pos- 
itive measurable function of cu. The discrete part Td{^^) which is a sum of 
Dirac functions, and the singular continuous part Tsd'-^) is simply the rest. 
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These different components are directly related to the corresponding compo- 
nents of the spectrum of the phonon bath. We will discuss this relationship 
only qualitatively, without attempting any rigor. The decomposition of the 
positive measure T{u!)du! into these three parts implies a corresponding de- 
composition T{t) — Tac{t) + ^d{t) + T sc{t) .Their qualitative time dependence 
will be discussed below for the absolutely continuous and the discrete part. 

(i) Absolutely Continuous Phonon spectrum 

For the eigenvalues cui, e Sac which belong to the absolutely continuous 
part of the phonon spectrum, the eigenvectors u^'^'' corresponding to cUi, are 
extended, their components u'f^ of the normalized state are proportional to 
1 / \/N where is the number of sites i of the harmonic part of our system 
which is supposed to go to infinity. Since we assumed that the coupling with 
the phonon bath Cj is local that is (12) is fulfilled, we can write 




where c,y is a finite number independent on N for N ^ oo. Moreover, 
the extended eigenvector u^^^ and consequently c^, depend on a continuous 

parameter (for example the wave vector). An important consequence is that 
the average of Cj, for uj,^ E [uj — 6uj, uo + 5uo\ in the limit of an infinite system 
and small 5u} is a well-defined measurable function of ui 

(c.).,.=a. = c(a;) = (c«(a;)) . (29) 

Thus, the contribution of the absolutely continuous phonon spectrum to the 
Fourier transform eq.(26) of the Kernel T{ijj)duj is nothing but its absolutely 
continuous part Tac{'^)- If A*(^) is the density of states of the absolutely 
continuous phonon spectrum the absolutely continuous measure is given by 

^aj{^)d^ = T l^{^)dt^ (30) 



which implies 

rSW = rS'(-t) = / , ^=(^cosM)M^)<i.. (3i) 

Since //(a;) is absolutely integrable, Tac{t) decays as 1/f or faster, due to the 
Riemann-Lebesgue lemma. This holds e.g. for a spatially periodic system in 
any dimension and for a three- or higher-dimensional system with disorder 
which is not too strong. 
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(ii) Discrete Phonon Spectrum 



For the eigenvalues cUj, e which belong to the discrete part of the phonon 
bath, u^''^ is a localized square summable eigenmode. Then, the scalar prod- 
uct Ci, — u^'^^^C is finite. Considering an arbitrary closed interval Xd C 5^, 
we require a stronger convergence than (14) for the discrete case: 

E ^ = < +00 . (32) 



Actually, this assumption should be fulfilled in most realistic physical situ- 
ations with a discrete spectrum. 

Generally, the eigenmodes w*^*^^ decay exponentially with respect to \i\ and 
there is a uniform bound on the localization length provided ujiy G Id (pro- 
vided this interval does not contain the mobility edges). If we assume the 
couphng constants Ci fulfill (12), then |C,^| = |w*^^)*C| decays exponentially 
as a function of di,, the distance from the "center of mass" of mode u. Then, 
since o;^ is bounded from below by a positive nonvanishing number, series 
(32) converges. Note, however, that even for a system with acoustic phonons 
the series (32) may exist because the translational invariance implies at least 
Ci, ~ cUi, for cui, — > 0. However, in principle, there is no need for exponential 
localization, the eigenmodes could decay with some algebraic law providing 
they are square summable. 

The contribution of the discrete spectrum of the phonon bath to the Fourier 
transform measure (26) is nothing but the discrete part Tdi^-^) which is a 
sum of Dirac functions 



This immediately leads to 

rSW=rS(-t)= E ^^cosa;.^ . (34) 

Because of (14), Td{t) is an absolutely convergent discrete sum of cosine 
functions. Consequently, in general it is an almost periodic function which 

does not decay to zero at infinity but oscillates between a lower and an up- 
per bound. The situation of a pure discrete spectrum occurs in disordered 
systems, in \d or 2d or in higher dimensions with strong enough disorder. 
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(iii) Singular Continuous Phonon Spectrum 

By definition, the singular continuous spectrum is the rest of the phonon 
spectrum after removing the discrete and the absolutely continuous part. 
If there is a singular continuous part in the phonon spectrum, it simply 
appears as the singular continuous component Tsc{^) of r(a;). In standard 
models which are spatially periodic or random, there is no singular continu- 
ous spectrum but there are special models [27,28] which are neither spatially 
periodic nor random (for example quasiperiodic) and where the whole spec- 
trum is singular continuous. In such models, for uj^ e Sgci the eigenvectors 
u^'^^ are not square summable though are not uniformy extended as for ex- 
ample planewaves. They have a behavior at infinity which is not universal. 
For example, they might slowly decay with power laws, or they might not 
decay but with intermittent and sparse peaks with possibly fractal distribu- 
tion and zero density in average. We shall not much discuss this exceptional 
situation in this paper. 



If the disorder in "id systems is weak enough a mixed situation occurs with 
an absolutely continuous part in a central interval and a discrete one in the 
two intervals above and below the mobility edges. Then, r(a;) is the sum of 
a discrete part and an absolutely continuous component with supports which 
arc the corresponding frequency intervals. When there arc acoustic phonons 
in a random model, the spectrum is still fully discrete in one and likely two 
dimensions but with a localization length which diverges at zero frequency. 
In three dimensions, the low frequency phonon modes are extended. At larger 
frequency beyond a mobility edge, the spectrum is discrete. In all these cases, 
the density of states remains a smooth measure indistinctly the spectrum is 
absolutely continuous or discrete. 

There are other situations we shall not consider here where the spectrum may 
be absolutely continuous but where the support is a fat Cantor set [29] (then 
r(ac)(<^) is not a continuous function of uj) or when the spectrum is discrete 
but is not dense on finite intervals (for example in smooth quasiperiodic poten- 
tials) or when the spectrum is purely singular continuous (quasiperiodic square 
potentials [27,28]). We shall not discuss here these exceptional situations. 

The Markov approximation is the simplest approximation which consists in 
assuming that the spectrum of the phonon bath is only absolutely continuous 
with infinite support and the coupling constants are uniform. Thus T{uj) — 
f ac(<^) = const and Ta,p{t) — — '^la,i3S{t). Then, eq. (22) becomes 

+ K'//(qW) + A^q(i) - AC(t) = (35) 
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which contains a damping term proportional to the velocities and a Langevin 
force. The solution of this equation at zero temperature, (C(^) = 0), always 
converge to a static solution corresponding to the system at equilibrium. This 
result is essentially the consequence of the absolutely continuous phonon spec- 
trum which extends from frequency zero to infinity. 



4 Energy Dissipation and Asymptotic Solutions 

Starting from an initial condition where the phonon bath is at rest, i.e. £,^(0) = 
and e,y(0) = the energy provided by the anharmonic subsystem to this 
phonon bath may be viewed as energy dissipation for this subsystem. In the 
general case, the rate of energy dissipation is found to be related to the Fourier 
transform of the velocity- velocity correlation function of the trajectory instead 
of being proportional to its square velocity as in the original Langevin case. 
At zero temperature when the fluctuating random force in eq.(19) vanishes, 
the possible asymptotic solutions are those which could persist without any 
energy dissipation. Since the dissipation rate depends on a time average, the 
role of the transient dynamics toward the asymptotic solution disappears. 

Actually, the asymptotic solution might not be always a static solution when 
the structure of the phonon spectrum is not uniform. We already know situ- 
ations when the phonon spectrum is absolutely continuous but bounded with 
an upper bound (and possibly with frequency gaps where the density of states 
vanishes, where (interband) DBs are asymptotic solutions. Then their funda- 
mental frequency (and harmonics) necessarily belong to gaps of the phonon 
spectrum [5,7]. 

We also know that when the phonon spectrum is purely discrete that there are 
time-periodic solutions (called intraband DBs) with fundamental frequencies 
which belong to a certain fat Cantor set (with finite measure) [19] dense in 
the phonon spectrum but which does not contain any phonon eigenfrequencies. 
However, we do not know if they could be asymptotic solutions for other initial 
conditions which are not those of that solutions. We would hke to comment in 
that section on the relationship between vanishing energy dissipation and the 
possible asymptotic solutions for arbitrary initial conditions. This relationship 
depends sensitively on the properties of the phonon spectrum. 

Let us first investigate the energy transferred after a long time to the phonon 
bath of a solution q{t) of eq.(19) (or of the equivalent equation (22)) at zero 
temperature, i.e. for ^(t) = 0. The consequence for the asymptotic solutions 
will be discussed in the final subsection. 
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4-1 Derivation of an Energy Dissipation Functional 



For time T > given, let us first define functions 



9a;T{i) = XT{t)qait) (36) 

where Xt(^) = f*^^ < t < T and Xt(^) = for t < and t > T, i.e. 

Xrit) is tlie normalized characteristic function on [0,T]. Its Fourier transform 
is 



r+oo 1 pT 

~9a;T{co) = / ga;Tit)e^'dt = -= / g„(i)e-*di (37) 

J—oo y 1 JO 

Let us now consider the effective energy of the anharmonic part of the system 
defined as 

E{t)^^mr + Veff{q{t)). (38) 

Using cq.(22) with <^{t) — 0, (since zero temperature is assumed) the time 
derivative of this energy becomes 

E(t) = -X^ f qit) ■ Tit - T)q(T)dT . (39) 

^0 

The average of this time derivative from time to T is 



= -iA'l r dt r g(t) ■ T(t - T)q(T)dT 

2 I Jo Jo 

1 r+oo r+oo 

^--X' dt grit) -Tit -T)gT{T)dT 

Z J—oo J—oo 



(40) 



where we used Va^pit — r) = r^^Q(r — t) and definition (36). Then, using the in- 
verse Fourier transform g^-rit) = 1/2% e-''^*ga-T{'^)du and e'^'^+'^'^^dt = 
27rS{u! + a;'), we obtain the energy dissipation functional 
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47r J-oo 



Definition (24) of T{t) implies that its Fourier transform f (a;) is a positive 

semidefinite matrix for any cu, Consequently we have ET{t) < 0. This result 

is physically obvious because the energy variation E{T) — E{Q) = TEt of 
the anharmonic part of the system has to be negative at any time T since 
the total energy must be conserved and there was initially no energy in the 
phonon bath. 

E{T) is obviously bounded from below (for physical stability of the model) and 
from above by the total finite energy which is initially injected to the system. 
Consequently, we have limT^+oo t ~ ^- Thus a necessary condition for 
the asymptotic solutions is that 



lim Et^O, . (42) 
r-»+oo 

The limit of the products in (41) for T — +oo might not be a smooth 
function of oo anymore. For example if q{t) is time periodic, it should ex- 
hibit Dirac peaks. Considering the limit within the topology of weak con- 
vergence of measures, this limit can be defined as a positive measure ma- 
trix Kf}Aijj)duj. Then, it is proven in appendix A, that as a measure, this 
limit is the Fourier transform of the velocity-velocity correlation functions 
Ka,p{r) = <ia{t)qf3{t + r) = limT-.+oo lo qa{t)qa{t + r)dt: 



Ka,p{u;)du; = Jim g*.j,{u;)gp.T(u;)du; = ( [ Ka,/3{r)e"^''dT) du (43) 

r^+oo \J-oo / 



4-2 Energy Dissipation for a Sinusoidal Periodic Driving Force 



For a given solution q{t) of eq. (22) with ({t) = we have derived (in subsec- 
tion 4.1) an exact formula for the averaged energy dissipation rate. We will 
apply this formula to a sinusoidal trajectory q{t) = acos{Qt — 9) with fre- 
quency Q and arbitrary phase 6. We have q(t) = —aQsm{Qt — 6) which can 
also be interpreted as a driving force. The normalized characteristic function 
fulfills / XT{t)dt = 1. Its Fourier transform 



22 



fulfills J \xTi^)\'^duj = 27r, liiiiT^+oo |xt(^^)P = 2-k5{lj) (for a; = 0, we set 
Xt(0) — VT) and lim XT(<^)Xr('^') = for a; 7^ a;. Then, we can write 

T— »oo 

eq.(37) as ga-^ri^) = ^ J Qai^J - uj')xT{oj')duj' where q{u;) = f q{t)e''^*dt = 
ivran (e-^^(5(u; + - e^^(5(u; - Q)). This yields: 

^a;r(c^) = (e-^'xT(c^ + - e^^XT(a; - Q)) (45) 

Then the average frequency dependent damping over time T for that time 
periodic function q{t) follows from (41) 



- (e"'''(XT(a; + Q)xt(^ - + cc. 



duj . (46) 



The average dissipation rate depends sensitively on the damping kernel 
r(a;) which itself depends on the phonon spectrum. We will discuss this de- 
pendence of Et for T — > 00 for an absolutely continuous and a discrete phonon 
spectrum. 

(i) Absolutely Continuous Phonon Spectrum. 

As discussed in (i) of subsection 3.2, T{u;)du; = Tac{u;)du; is an absolutely 
continuous measure where Fad^^) is a positive continuous function on its 
support which is the absolutely continuous spectrum Sac- It has been as- 
sumed to be the union of one or several intervals. Tac{'^)duj is vanishing 
elsewhere. Then, using the properties of Xt('-^) from above, we get a result 
independent of the phase 9 



-C/no — — 



j «af a,/3(<^)a/3 {5{uj + Vt) + 5{uj - Q)) dhJ 

a,f3 

= -(^)'«*r(Q)a (47) 

where i'i—u;) — T{u;) has been used. 

Consequently, there is energy dissipation, only, if the driving frequency fl or 
(— f2) belongs to the absolutely continuous spectrum Sac- In that case the 
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rate of the energy dissipation is well-defined and the energy of the phonon 
bath linearly increases in average with respect to time. 

If Q belongs to the edge of one of the intervals the Dirac delta functions in 
eq.(47) will contribute only half to the integral, i.e. is still a constant. 
Accordingly, £^7- ~ T for T — > oo. However, this is not true anymore if 
in one-dimensional systems there is a van Hove singularity. In that case 
T{uj) diverges and a superdissipative behavior should be expected(see next 
subsection) . 

When f2 ^ Sac we have min^^g^,^^ |f2 — a;| = daciP) > and consequently 

Ichiriuj ±Q)\ < Then it comes out from eq.(46) that \Et\ < 

K/T where K is some finite constant. Since Et = {E{T) — E{0))/T, the 
effective energy E{T) is monotonously decaying towards hm^^+oo E(T) = 
£■00, but remains bounded from below at infinite time. 

(ii) Discrete Phonon Spectrum 

This case has been discussed in (n) of subsection 3.2. Ta{(^)duj is a discrete 
measure (33). Substituting this into eq.(41) and taking eq.(45) into account 
yields 



Et = — J^a^le-^^XTK + ^) - e^'xrK - (48) 

where Ui, = a^CuCla/ul = ciaCu,aCu,i3(ii3/'-^l are positive coefficients. 

a,f3 

Since the positive series J2u (^u < +00 is absolutely convergent, we can label 
the countable set of eigenfrequencies by positive integers u & J\f such that 
the sequence ai, is monotone decaying (oi, > aj^+i for any u ). Thus, the res- 
onances u are simply ordered according to their strength. Then, a^, can be 
supposed to converge exponentially to zero as in standard physical models 
with discrete spectrum. The f2-dependent energy dissipation will be more 
subtle, compared to case (i) of an absolutely continuous phonon spectrum. 

When fl ^ Sd does not belong to the discrete spectrum (which is the clo- 
sure of the set of eigenvalues), the same proof as above for Q ^ Sac yields 
again that there is no energy dissipation and that the phonon bath energy 
remains bounded from below. When Q e iS^ and when there is a resonance 
that is there exists uq such that we have precisely Q = ±<^i/o> then series 
(48) contains a term Xt{^ i ^vo) — VT which diverges when T — -|-oo. 
Consequently, E^ diverges as —T when T — ^ +00 while the energy absorbed 
by the phonon bath E^ diverges proportionally to T^. This result should be 
expected since we know that the energy of an harmonic oscillator driven at 
its resonant frequency diverge as a function of time as t^. 
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However when Q & Sd but Q ^ ±uj,j Vz/, the situation is more subtle because 
series (48) though it does not contain any diverging term, contains terms 
with denominators | {uji, ± fl) \ which may be arbitrarily small and may cause 
divergency. Inspired by the construction of the Liouville numbers from ra- 
tional numbers, we construct a subset of Sd which has full measure and 
where this series is absolutely convergent for all fl in that subset. 

With Xr(<^) from eq.(44) we get from eq.(48) 

F(T) = E{T) - E{0)=TEt = X 
uj,y + Q, uj,^ — fl 



F{T) is an infinite series of periodic functions of T. If this series is ab- 
solutely convergent, its sum is an almost periodic function of T (in the 
sense of H.Bohr [30]) which is smooth and bounded. A sufficient condition 
for absolute convergence of series (49) is that the positive series A{fl) — 

( KW + Ji^^:^) convergent. 

In appendix B we will prove that there exists a set C of zero Lebesgue 
measure such that the series A is convergent for all fl e TZ\C. This imphes 
that there exists a set TZ\C of full measure such that the energy —F{T) 
transferred to the phonon bath remains finite and bounded for all times, 
provided Q belongs to TZ\C. On the other side, it is also proven in appendix 
B that there exists a subset C of £ of zero Lebesgues measure such that 
F{T) diverges for T — > oo, if Q e £,'. In addition it is proven that under 
some conditions for {uiy} an arbitrary eigenfrequency belongs to C with 
probability zero, only. This excludes F{T) ~ for all fl G C. For fl G 
TZ\C it is F{T) ~ and F{T) ~ T for an ordinary dissipative Langevin 
bath with an absolutely continuous spectrum. Therefore we may conjecture 
that F{T) ~ if Q G with an exponent a which may depend on fl. 
Compared to ordinary dissipation where a — 1, situations where 1 < a < 2 
may be called superdissipative and those where < a < 1 subdissipative. 

It has been shown that the classical dissipation kernel T{t) determines the 
nonlocal action occurring in the path integral representation of the quan- 
tum mechanical propagator, after elimination of the harmonic d.o.f. [31]. 
The nonlocal action describes quantum dissipation. Depending on the low 
frequency behavior of the spectral density J{uj) (corresponding to r{uj)), 
i.e. J{uj) ~ uj'"^, for uj ^ 0, qualitative different quantum behavior is found. 
For instance, for a particle in a double well coupled to a harmonic bath 
there can be ohmic damping (s = 1), superohmic (s > 1) and subohmic 
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(0 < s < 1) damping [32]. Since the Q-dependent dissipation discussed 
above is not related to the low frequency behavior of T{uj), the implication 
of the subdissipative and superdisspative behavior for quantum dissipation 
is not obvious. 

In conclusion, we have obtained here a rigorous proof that when the 
phonon spectrum is discrete and with the assumption of fast enough decay 
of coefficients (which is generally fulfilled in physical models), there is 
no energy dissipation by the phonon bath when it is driven (locally) by 
a time periodic force for most choices of its frequency Q (i.e. probability 
1), when Q belongs to the phonon spectrum Sd (the closure of the set of 
eigenvalues). This result first sharply contrasts the properties of a Langevin 
bath with absolutely continuous spectrum which is always dissipative and 
second, gives evidence for the existence of IDB. We also proved the existence 
of a zero measure subset C of frequencies for the driving force where the 
energy absorbed by the phonon bath is unbounded as a function of time. In 
that case, the phonon bath may be superdissipative or subdissipative. 

(iii) Mixed Case 

When the phonon spectrum contains both an absolutely continuous com- 
ponent and a discrete component, energy dissipation is simply given by the 
sum of these two contributions. Thus for having no energy dissipation, Q 
has not to belong to Sac but may belong to the subset 7l\C in Sd defined 
above where there is no energy dissipation. Of course, there is no energy 
dissipation when Q does not belong to any components of the spectrum 
Sac U Sd- 

4-3 Possible Asymptotic Solutions 

If one chooses a general initial condition q(0), q(0), {e^^O)}, {^(O)} for the full 
system one may ask for the asymptotic behavior of the solution x{t) = (q(t), 
{£u{t)}) of eq.(22) for i — > oo. The asymptotic solutions of a trajectory 
X(t) = {q{t), {e^{t)}) generated from an arbitrary initial condition are de- 
fined as the trajectories in the intersection PIt ^ for all time T of the closure 
of the trajectory Tt = X{t)t>T starting from time T. This set could be a 
fixed point corresponding to a static solution. It could be an invariant circle. 
Then the asymptotic trajectory would be time periodic (Discrete Breathers). 
It could be also an invariant torus and then the asymptotic trajectory would be 
quasiperiodic. We should also not exclude a priori situations were the limit set 
is a strange attractor with chaotic trajectories. We wish to discriminate about 
the possible limits which are allowed in principle according to the properties 
of the phonon bath. 
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In the following we will discuss qualitatively which asymptotic solutions may 
be possible. From the average energy dissipation rate (cq.(41)) we get a nonover- 
lapping criterion for the support of both positive measures r{u!)du! and K{uj)duj. 
This criterion, together with the generation of higher harmonics, can be used 
to specify features of the asymptotic solutions. 

(i) Nonoverlapping criteria 

Equation (41) relates the rate of energy dissipation after a finite time T 
to the smooth function gri^jj)- For T — > oo, gri^^) is related to the Fourier 
transform of the velocity-velocity correlation function (43) which does not 
depend on the possible transient regime near the initial time but essentially 
on the asymptotic solution. For T ^ oo, the rate of energy dissipation be- 
comes 

— — , r+oo _ 

= hm Et^-—Y. r^,f,{u;)Kf,,^{u;)dw (50) 

T-.+0O 47r J-oo 

The positive matrix of limit measures K{uo)duo = dp{uo) (Fourier transform 
of the matrix of velocity- velocity correlation functions) can be decomposed 
as usual in three components (absolutely continuous, discrete and singular 
continuous) 

dp{ix}) = dpac{io) + dpd{oj) + dpsc{oj) 

Then, strictly speaking, the integral in eq.(50) of the product of two pos- 
itive generalized functions may be mathematically undefined when both 
measures T{u})du} and K{u})du} are not purely absolutely continuous. Such 
a situation may occur for example when both of these measures exhibits 
a Dirac component at the same frequency lOv However, in that case, the 
physical situation is clear because the anharmonic system would drive the 
phonon bath with a resonant time periodic force component at frequency 
ujy. Then, the energy which is provided to the resonant harmonic oscillator 
should grow as a function of time as t^, faster than t. This resonant situation 
would correspond to a superdissipative behavior where l^'ool — +cio since 
the energy dissipation per unit time would diverges as t ioi t ^ +oo. More 
generally, we have to ascribe an infinite positive value to the integral (50) 
if its integrand is a product of two singular continuous or discrete functions 
and their supports overlap. 

Since the support of a positive measure plays a crucial role in the fol- 
lowing and not all readers may be familiar with it we give a brief defi- 
nition. A support /Cac for the absolutely continuous part dpac{oj) can be 
chosen as the set of points where the derivative is well-defined (for 

most uj) and is nonvanishing (it is positive). It is a measurable set with 
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nonzero measure. A support Kd (the smallest) of the discrete component 
5pd{uj) = J2n kn^iyj — VLn)du is the countable set of points where /c„ > 0. 
The support fCsc of the singular continuous measure dpsd^) can be chosen 
as the set of points where '^^^l^^ diverges to infinity (This limit is defined for 
most u and is zero by definition). It is a measurable set which is uncountable. 

We assume that the support of Ta^/s{uj)duj is the same for all Eq.(26) 
shows that this holds if for given u the coupling constants C,y^a are nonzero 
for all a. This may be not true in presence of special spatial symmetries. 
The same assumption is made the for the support of Ka^p{uj)duj. In generic 
models all variables qa{t), « = 1, 2, 3, ■ ■ ■ are coupled one with each other. 
Consequently, we assume that the support of the Fourier transform of the 
correlation functions qaif)<ii3ii + t) does not depend on This implies 
that {Kai3{uj)d{ijj)} have a common support. 

The main result is that since the initial energy of the anharmonic system is 
finite and the total energy is bounded from below (for stability reasons), a 
necessary condition for the asymptotic solutions, is that the average rate of 
energy dissipation at time infinity given by the integral (50) be strictly zero. 
This implies the support of the positive measures T{u})dui and K{ui)duj do 
not overlap. 

(ii) Conjectures for the limit solutions 

Since gri^^) for T — ^> oo is not the Fourier transform of the velocity q{t) 
which itself is directly related to the Fourier transform q{uj) of the trajec- 
tory q{t), the support of K{uj)duj may not yield information on q{uj)duj. 
If q{t) is periodic, quasiperiodic or almost period [ll then K{t) is periodic, 
quasiperiodic or almost periodic. In this case the support of K{u)du! is 
identical with that of q{uj)duj. However, if q{t) is chaotic it is not clear how 
q{(jj)duj is related to K{uj)duj. 

If VL is in the support of q{uj)duj, the nonlinearities in eq.(22) will generate 
all higher harmonics with frequencies Qp, with p an integer. Similarly, if Qi 
and Q2 are in the support, the nonlinearities will generate higher harmonics 
with frequencies niQi + 77,2^^2; with ni and n2 integers. We assume that this 
also holds for the support of K{uj)duj. Therefore we conjecture that because 
of nonlinearities of the anharmonic system, the supports Kac (resp. /C^, /C^c) 
have the following invariance properties. If Vti G /Cac and 0.2 G /Cqc (resp. 
A^d, ^scj) then we have riiVti + 712^2 G ICac (resp. /C^, Ksc) for any pair of 



Almost periodic functions are convergent series of sine and cosine functions with 
a countable set of frequencies. Quasiperiodic functions are a special case of almost 
periodic functions where the involved frequencies are linear combinations of a finite 
number of frequencies with integer coefficients 
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integers ni and n2- A measurable set which has this property is either the 
full real axis or has a zero measure (proof in appendix D). Consequently, 
JCac should be the whole real axis and therefore should overlap the support 
of the phonon spectrum S. This generates energy dissipation. Consequently, 
there should not exist any absolutely continuous component in the Fourier 
spectrum of the velocity-velocity correlation function but only a discrete 
component or a singular continuous one (or both). 

In case of a discrete component there are two possibilities. First, K{uj)duj ~ 
5{uj — Q)dijj, i.e. K{t) is periodic with period T — 27r/Q which is consistent 
with an asymptotically periodic solution q{t + T) = q(t) for t oo. Second, 
K(uj)duj = Y.kn^ioJ — fln)du! with X^^n^n 7^ for all integers {m„}, i.e. 

n n 

K(t) is quasiperiodic or almost periodic. This again is consistent with an 
asymptotically quasiperiodic or almost periodic solution q{t). 

In order that an asymptotically periodic solution exists its frequency has 
not been in the support of r{uj)du. This requires f2 7^ u;^ for all u. We 
have already shown (cf. (ii) of subsection 4.2) that this condition is not 
sufficient. This point will be discussed in great detail in section 5. For an 
asymptotically quasiperiodic or almost periodic solution a necessary condi- 
tion is fi„ 7^ uji, for and n and all i/. We do not discuss this case but mention 
that the existence of such solutions were proven for a special class of models 
with quenched disorder [22]. We expect that such solutions also exist for 
the present model. 

Since solutions which are neither periodic, nor quasiperiodic (or almost 
periodic) are chaotic, the case of a singular continuous component of K{uj)djj 
can only be related to asymptotically chaotic solutions. The fact that the 
component is singular and not absolutely continuous (which has already 
been excluded) may imply that the chaos is not fully developed. This subtle 
point is also not discussed further. 

Finally, the asymptotic solution can also be a constant one, i.e. lim qit) — 

t—fOO 

q*, where q* is a local minimum of Vefi{q). In that case it is qoo{uj)duj — 
27rq*5{uj) and K{ijj)dw ~ uj'^5{uj)duj, s.t. E^o — I i3{i^)Ki3 a{uj)du; — 0. 

Here we used that Ta/siO) exists. 

This discussion shows that asymptotic solutions can be constant, periodic, 
quasiperiodic (or almost periodic) and may be weakly chaotic. As far as we 
know the existence of asymptotically weak chaotic solutions of eq.(22) has 
never been investigated. We should also stress again that we have discussed 
possible asymptotic solutions, only. Wc have no information which of these 
asymptotic solutions will occur for given initial conditions, e.g. q{0), q{0), 
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{£,(0) = 0}, {£,(0) = 0}. 



5 Existence Proof of Intraband Discrete Breathers 

In the previous section we have argued that the solution q{t) of eq.(19) (or 
equivalently eq.(22)) may converge towards an asymptotic solution which 
might be periodic, quasiperiodic or weakly chaotic. 

We prove now that indeed there may exist exact periodic solutions under some 
conditions. When the phonon spectrum has an absolutely continuous compo- 
nent, there may exist exact time periodic solutions providing their frequency 
and higher harmonics are outside the absolutely continuous spectrum <Soc- 
When there is a discrete component, those solutions may exist with frequency 
and higher harmonics inside the discrete spectrum Sa but sufficiently different 
from the phonon frequencies cuu- In the present section, we prove the existence 
of such periodic solutions. 

Considering a time periodic solution of the uncoupled (A = 0) Hamiltonian 
7fo(p, q) with frequency fl = 2n/T and period T, we ask whether this solution 
persists as an exact time periodic solution at the same frequency when the 
coupling A is turned on at least up to a nonvanishing value, will be considered 
as a given number which has to be chosen in an appropriate set in order this 
solution persists. 

We also need to assume some conditions concerning the time periodic solutions 
which are generally fulfilled in order to answer this question. This time periodic 
solutions should belong to a one parameter family of time periodic solutions. 
Thus, this family can be represented locally in action angle variables. Its energy 
Hb{I) only depends on its action / (at least defined in some interval), then 
i?^(/) is the frequency of the solution. We need to assume that this solution is 
not marginal at the considered frequency fl that is H"{I) ^ when H[j{I) = Q. 
The possible existence of such time periodic solutions in complex systems is 
already known and for example may be based on earlier proofs of (interband) 
DBs (for which the fundamental frequency and its harmonics are not resonant 
with the phonon band) considered in [10]. 

The equation of motion for a periodic solution q(t; Q) with period T = 27t/Q 
of the anharmonic subsystem coupled to the collection of harmonic oscillators 
is (see appendix D) 

T 

q(t; n) + 1/'(q(t; O)) - J G{t - r; f))q(T; n)dT = (51) 
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with 

GaAf, ^) - E o cosa;.(t - 0<i<T (52) 

V 2a;,,sm(a;,,-|-) 2 

which should no be confused with Ga,p{t). Gap{t; fi) is periodically continued 
for i < and t > T. A necessary condition for the existence of IDB with 
frequency fl is that the periodic kernel Ga,i3{t; ^l) is a bounded and smooth 
function in t for given fl. This condition will be discussed in the following. 



5.1 Convergence of Kernel Contribution due to the Absolutely Continuous 
Spectrum 



According to the hypothesis, eqs. (28), (29), the absolutely continuous part 
becomes the integral of a smooth function of cu in the limit of N to infinity: 



Mac) 



{t,n) 



COS OJi, 



Ca{uj)Cf3{u) 
2U! SiUTT^ 



T 

T 



COS a; 



(t - n{uj) duj 



(53) 



which diverges for cu — pQ, if pfl e Sac for some nonzero integer p. Then, 
G^"''^\t; Q) is not a smooth bounded function of time. Note, the divergence 
at a; = (p = 0) is harmless because of the hypothesis eq.(14) of finite 
reorganization energy. Thus, the contribution of the absolutely continuous 
spectrum, G^"''^\t; Q), is a smooth bounded function of time when 

Q^U% ■ (54) 

p>0 P 

This is the standard condition expected from earlier theories for the existence 
of DBs which assume the absence of resonance of the fundamental frequency 

and its higher harmonics with the absolutely (continuous) phonon spectrum. 
It is consistent with the conclusions drawn from eq.(46) in (i) of subsection 
4.2. 



5.2 Convergence of the Kernel Contribution due to the Discrete Spectrum 



The contribution due to the discrete spectrum 
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cannot be transformed into the integral of a smooth function. Eq.(55) clearly 
demonstrates the problem of small denominators which makes the proof of 
the existence of IDB with frequency Q in presence of a dense discrete phonon 
spectrum nontrivial. This problem of small denominator is quite similar to 
that we encountered in (ii) of subsection 4.2. In the following we will study 
this small denominator problem. 

For given constant c > and given Q, we define a closed interval or a finite 
union of closed intervals Id strictly included in the discrete spectrum Sd such 
that for uj G Td, \ sin vr^l < c. c is arbitrary except that it must be chosen small 
enough in order Ta does not contain the edges of Sd (mobility edges). This set 
Id contains the values of e Sd- Then, for cui, e Sd\Id, the denominators 
I sinvr^l > c > are bounded from below by c. 

We split again the discrete kernel (55) into two parts G'^'^^(t; = G^'''^\t; + 
G^'^'\t;n) where 



2a;^,sin7r^ V 2 

is a smooth bounded function of time because of the absence of small de- 
nominators (they are always larger than cui, and therefore the sum converges 
because of the hypothesis of finite reorganization energy (14); the divergence 
at o;,^ = is harmless since Ci,^a '-^ cui, for cu^ ^ 0, as already mentioned before) 

The small denominators are in the second part 



In order to insure the convergence of (57), we choose f2 in a particular subset 
C{k) of the real positive axis defined such that for all uji, e Sd, we have 



< A;i^ . (58) 





|C.| 


2 




sinvr^l 



Then, if Q e C(/c), we have 



32 







2 







which is convergent according to the hypothesis (32). Consequently, for Q e 
C{k\ the discrete spectrum contribution Qi'^'^{t\ Q) to the kernel G(i; Q) in eq. 
(55) is a smooth bounded function of time. 

We now have to prove that C(A;) is not an empty set if we choose the parameter 
k in an appropriate way. First, we choose k large enough in order < e,^ < | 
can be defined (uniquely) for all lo^ e by 



2k 

Then, condition (58) becomes 



'QJ\= sinvrei, . (60) 



which requires for all p > and all lOv G Sa 



siuTrei, < |sin7r— ^| (61) 



O ^ V{k) = u 



(62) 



C{k) is the complementary set C{k) = TZ^\D{k). The total measure of the 
gaps of C{k) can be easily bounded from definition (62) for u^, e Id where Id 
is any interval or union of intervals strictly included in Sd as 

^mk)) < E E 4^ < E 2e.a;e f E -AtiI ■ (63) 

i^^eTdP>oP ^i' ui^eSd \p>qP ^l^j 



Since < e,, < |, we have < 2ei, < siuTrej, = ^|Ci/| < , using 

definition (60) and inequality (58) . Then using condition (32), we obtain 



l,{V{k)) < jB{Id) . (64) 



where 



a;3 



Choosing k large, the Lebesgue measure of T>{k) becomes small so that the 
measure of the Cantor set C{k) tends to full measure in Up>o ^• 
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5.3 Proof of existence of IDB 



After having proved the existence of a set C{k) s.t. Ga,p{t]VL) is a bounded 
and smooth function of e C{k) we now apply the imphcit function theorem 
to prove the existence of IDB. 

We consider the Banach space Ci of T-periodic different iable functions q{t)) 
with continuous time derivatives which are time periodic p{t+T) = pit), q{t + 
T) — q{t) and fulfill time reversibility: p{t) — —p{—t) and q{t) — q{—t). The 

norm in jCi is defined as ||p(t), qf(i)||i = supJ |p(i)|, \qit)\, 



\p{t)\,\q{t)\j. We consider the Banach space Cq of T-periodic continuous 
and time reversible functions {v{t),w{t)) with the norm 

s^Pt\\'"i't)\AMt)\\- 



Let us consider for a given Q G C{k), the map T{p, q; A) defined as 



V 



w 



'r{p,q;X) 



yP + V\q) - X' /o' G{t - r; l])q(r)dr j 



(65) 



It maps Ci in Cq and is continuous and different iable with respect to {p, q) 
and A. The zeros of T(p, qr; A) yield time periodic solutions of eq. (51). 

At A = 0, let us assume that the anharmonic subsystem has a periodic 
and time reversible solution qo{t) with period T — 27r/Q with Q, e C{k). 
If dT{p, q; 0) is invcrtiblc at {po, qo) then this solution can be continued in A 
in a neighborhood of zero up to a nonzero value. We have 



dT{po,qo;0) 



^ Sp\ ( 5p — 6q 
^Sq) [5p+V"{qo)5q^ 



(66) 



The condition for invcrtibality of dT{po, qo; 0) follows from the Floquet anal- 
ysis [7]. Setting the r.h.s. of (66) to zero yields the linearized equation of 
motion with respect to the periodic and time reversible solution {po{t), qo(^)) 
of the anharmonic subsystem. Integration with respect to time over one pe- 
riod leads to the Floquet matrix of the unperturbed DB. Then, dT{po, qo] 0) 
is invertible, if the Floquet matrix has only a single eigenvalue which equals 
one. Since its corresponding eigenfunction is antisymmetric with respect to 
time reversion, this eigenmode is harmless. 
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6 Summciry and conclusions 



We have assumed that a general macroscopic classical particle system with 
arbitrary interactions can be decomposed into a microscopic or macroscopic 
anharmonic subsystem with coordinates, q = (g^), and anharmonic interac- 
tion energy V{q), and a remaining macroscopic part with harmonic interac- 
tions. In order that the linear spectrum may possess a continuous part one has 
to take the thermodynamic limit of the harmonic subsystem. The interaction 
between the harmonic displacement variables Ui at site i with is assumed 
to be bilinear with coupling constants AC^.q. A is introduced for technical con- 
venience. A realistic system for which these assumptions may apply is, e.g. 
a complex impurity embedded in a lattice where the latter can be treated in 
harmonic approximation. The harmonic part can be diagonalized by introduc- 
ing normal mode coordinates u^, for eigenmodes with frequency uj^,. Then our 
model takes the form which we call a Zwanzig-Kac model [3,4] for which the 
harmonic degrees of freedom can be eliminated. This results in the integro- 
differential equation eq. (19) for q{t), which is of Langevin type with memory 
kernels G(t). The properties of G{t), or equivalently of the dissipation kernels 
T{t) — —G{t), depend strongly on the spectrum of the harmonic part. 

The elimination of a macroscopic number of degrees of freedom leads to fluc- 
tuating forces and (as always believed) to dissipation. The former vanish at 
zero temperature. Then it is the common belief that q{t) for arbitrary initial 
conditions (g(0), q{0)) converges to one of the local minima of the effective 
potential Vcs{q)- However, this is only true if the phonon spectrum of the 
harmonic part is absolutely continuous with infinite support. In general, its 
support may be not infinite and in addition there may also exist a discrete 
spectrum, e.g. due to disorder. It has been one of the major goals of our 
contribution to prove that energy dissipation strongly depends on both, the 
specific properties of the phonon spectrum and of the trajectories q{t). For 
rather general conditions we have succeeded to derive an exact formula for 
the energy dissipation rate Et of the anharmonic system averaged over a 
time interval of length T. For T — > oo it only involves Ta,i3{uj)Kp,a{'^): the 
product of the Fourier transform of the dissipation kernels ^a,/3{t) and of the 
velocity- velocity correlation functions Ka,p{t). For a finite initial energy and 

a Hamiltonian bounded from below it must be Et=oo = 0. This is possible, 
only, if either Ka,i3{uj) = 0, which holds for q{t) q^ (local minimum of 
Ves{q)) for t — >• oo, or the support of the measures T{u!)d{u!) and K{uj)dijj 
have no overlap. Besides the convergence of q{t) to one of the local minimum 
q* of Ves{q)i the nonoverlapping criteria allows possible asymptotic solutions 
which are periodic, quasiperiodic (or almost periodic) or may be even weakly 
chaotic. 

If weakly chaotic asymptotic trajectories really exist, the corresponding mea- 
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sure K{uj)duj must be singular continuous. An absolutely continuous measure 
part of K{uj)duj is ruled out, since due to a theorem (proved in appendix C) 
its support must be the full real line, which always has an overlap with the 
support of T{u)du. 

Quasiperiodic (or almost periodic) aymptotic solutions can not exist as well, 
if rac{i^)duj has an absolutely continuous component . Because, with Qi and 
in the support of K[uj)duj also riifli + n2^2 with Ui integers will be in 
that support. Then there exist an infinite number of integer pairs (ni, , 
such that niQi +^2^12 has an overlap with the support of Taci^^)du!, provided 
Qi/fl2 is an irrational number. Quasiperiodic (or almost periodic) solutions 
may exist if r{u)duj has a discrete support, only. We do not have a proof for 
this. For a special class of models with disorder the existence of such solutions 
were proven rigorously [22]. The model studied in the present paper, may 
posses almost periodic solutions, as well. 

Periodic asymptotic solutions can occur when rac{^^)duj has a compact abso- 
lutely continuous, a singular continuous and/or a discrete component. Their 

existence is supported by the calculation of the dissipation rate Et for a si- 
nusoidal force with frequency fl. In that case a zero measure set C exists (see 
appendix B), such that \Et\ is bounded for all e T^\^, the full measure 
complement of C. In that case a Fourier component of q{t) with frequency 
Q would be undamped. That dissipation might be even more subtle has been 
demonstrated by the proof of the existence of a subset C of C (see appendix 
B) . The zero measure set C does not contain any of the phonon frequencies cu^, 
with probability one (assuming uji, and uji,/, v ^ i/' as uncorrelated). Although 
Q & C is not in resonance with any phonon frequency there is dissipation. 
Depending on Q there may exist subdissipative and superdissipative behav- 
ior, i.e. Et increases, respectively, slower and faster than T, where the linear 
T-dependence is the behavior of ordinary dissipation. Whether Et follows a 
power law with < a < 1 (subdissipative) and 1 < a < 2 (superdissipa- 
tive), as conjectured in (ii) of subsection 4.2, is not obvious. Anyway, we have 
shown that dissipation can be rather complex. 

The existence of nonrelaxing solutions also shows that the phonon bath can 
act as an anomalous thermostat. Depending on the initial conditions for q{t), 
the perturbed anharmonic system may not relax to equilibrium, provided the 
initial energy is finite, i.e. of order in the total number of d.o.f. There may 
even exist a sharp transition between relaxing and nonrelaxing situations. As 
shown recently [33] for a reduced Fermi-Pasta-Ulam model with a single an- 
harmonic bond embedded in an infinite harmonic chain there exists a critical 
excitation amplitude Ac of the anharmonic bond such that the initial con- 
figuration relaxes to equilibrium for A < Ac and to a discrete breather for 
A > Ac- That such an incomplete energy spreading can occur has also been 
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proven for a general discrete nonlinear Sclirodinger equation in any spatial 
dimension and without and with disorder [34]. 

Another important goal of our paper has been to prove the existence of IDBs. 
The existence of IDBs with frequency Q has been shown by proving the 
smoothness and boundedness of the periodic kernel G{t; Q) with respect to 
time t. In that case, following reference [10], the implicit function theorem 
can be applied to the time periodic solutions at frequency fl of the uncoupled 
subsystem, i.e. at A = 0, in order to prove the existence of time periodic so- 
lutions at the same frequency Q for a nonzero value of A. The boundedness 
of G{t; fl) is not obvious due to the occurrence of small denominators (cf. 
eq. (52)). These small denominators make more explicit the resonance prob- 
lem between the frequency Q and the harmonic mode frequencies uUi, (which 
is also the basic problem, discussed in Ref. [15,16,17]). It is quite analogous 
to the small denominator problem in secular perturbation theory in classi- 
cal mechanics. We have proven that there exists a Cantor set C{k) of finite 
Lebesgue measure such that G{t; Q) is smooth and bounded for all Q G C{k). 
This implies that the periodic solutions of the isolated anharmonic subsystem 
with frequency fl e C{k) can be continued up to a small but finite interaction 
strength between the anharmonic and harmonic subsystem. Therefore, there 
arc IDBs for all fl e C{k), provided A is small enough. The approach in Ref. 
[15,16,17] leading to the same conclusion for a different class of models, is 
partly complementary to ours. The authors of Ref. [15,16,17] consider a fully 
anharmonic system with quartic anharmonicity whereas our model consists of 
an anharmonic subsystem embedded in a harmonic lattice. The anharmonic- 
ity can be rather general. For our class of model it is the elimination of the 
harmonic d.o.f. which finally leads to the small denominators in its explicit 
form of eq. (52). 

There are several open problems which could be studied but which go be- 
yond the present contribution. First, the linear stability of the IBDs is one 
of them. Second, the existence of quasiperiodic or almost periodic solutions, 
which have already been proven for a special model [22], and particularly of 
weakly chaotic solutions of the present model could be investigated. Third, 
there is the question whether an arbitrary initial condition (q(0), q{0)) with 
large enough energy does relax completely or not. If it does not, where does 
it converge to? Since almost periodic solutions (which correspond to the ex- 
istence of KAM tori) should be linearly stable, i.e. all Lyapunov exponents 
should be zero, one should exclude convergence to true tori. Thus, the limit 
solution could only be at the "edge" of the KAM tori region with still discrete 
or perhaps singular continuous spectrum? Fourth, if both scenarios, i.e. com- 
plete and incomplete spreading (depending on the initial condition) exist, it 
will be interesting to study the existence of a sharp transition between both. 
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A Relation between the time Fourier transform of a bounded func- 
tion and the Fourier transform of its autocorrelation function 



We prove that (if defined) the hmit of the square modulus of the Fourier 
transform (eq.(37)) of Qa^rit) (eq.(36)) is identical to the Fourier transform of 
the matrix of autocorrelation functions of the velocities, but in the measure 
sense only. Since these Fourier transforms might not be smooth functions, 
they should be considered as positive measures which may have a discrete or 
singular continuous part. In that case a direct analytical treatment would lead 
to convergence problems because the kernel of the involved integrals, may not 
vanish at infinity. It is thus worthwhile to present a more rigorous treatment 
showing a relation not in terms of smooth functions, but of measure. For 
simplicity, we consider here only the case with a single component for the 
velocity vector but the proof would readily extend to cases with arbitrarily 
many components. 

Let us consider a smooth bounded function q{t) (which could be for example 
the time dependant velocity we consider above) and its Fourier transform over 
a finite time T weighted by 1/ VT 



This is a smooth function of oj for any finite T. Its square modulus defines an 
absolutely continuous positive measure dariu!) — \gTii^)\'^du. We must now 
assume that limr^+oo (Tt{^^) = o'(^) is defined in the sense of weak convergence 
of these measures [25]. The limit is a priori any kind of measure absolutely 
continuous, singular continuous or discrete. This assumption means by defini- 
tion that for any bounded continuous function /(a;) (for the uniform topology) 
we have limr_^.+oo/ /(<^)(^crr(<^) = / f{'-^)d(T{uj). 

Let us consider the subset of bounded continuous functions f{uj) with a com- 
pact support for ensuring the existence of a smooth Fourier transform F{t) ~ 
J f{uj)e~'^*du which fulfills the condition / \F{t)\dt < -Foo and / \ f{uj)\du < 
-\-oo. Then, 



Since \q{t)\ is bounded, this double integral is absolutely convergent in the 
limit T — > -|-oo, it readily comes out with the new integration variables t and 




Iff 

T Jo Jo 




q{t)q{t')F{t-t')dtdt'. 
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T — t — t' that for any Li smooth function F{t), we have 



^lim^J fiu;)daTiu;) = / f{u;)da{u;) = F{t) (^^lim^ ^ q{t)q{t + T)dt ^ 

(A.l)' 

Consequently, for most values of r, the autocorrelation function 



C(r) = C(-r) = lim 1 T + r)di (A.2) 

r^+oo 1 Jo 

is a well-defined function (since |C(r)| is obviously bounded). Defining, the 
Fourier transform measure dc{u!) = (/ C{T)e^'^'^dT)du!, it comes out 

I f{uj)da{uj) = I f{uj)dc{uj) 

for all continuous bounded functions f{uj) with compact support. We can 
now take sequences of continuous function with compact support converging 
toward arbitrary continuous functions without compact support which proves 
the equality remains fulfilled for arbitrary continuous functions f{uj) Thus, 
we prove the equality of the Fourier transform of the autocorrelation function 
dc{u}) considered as a measure with the limit at infinite time of the modulus 
square of finite time Fourier transform of the trajectory averaged over time 
d(7{ijj). 

It appears from eq.A.l that assuming daT{uj) has a limit measure for T +00 
or that the autocorrelation function (A.2) is defined are equivalent conditions. 



B Derivation of set C and C 



A sufficient condition for the absolute convergence of series (49) is the conver- 
gence of 

Given an arbitrary convergent positive series {sy\ with Sj, > which decays 
slower than ai, for example as a power law of u such as s^, = (with e > 0) 
and given some positive integer A;, we can make that scric ^(^2) convergent if 
we choose Vt such that |r2 + a;i,| > r^/k and 1^2 — a;,^| > rj,/k with r,, = {j^Y^'^. 
This condition means ^ Ck where Ck is the union of open intervals 
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(B.2) 



The Lebesgue measure of this union of intervals is bounded as 

Since is supposed to converge exponentially to zero and since has been 
chosen to go to zero as a power law, series S = J2u^u = {j^Yl^ < cx) is abso- 
lutely convergent for chosen in the complementary set of Consequently, 
we have f^{Ck) < 4:3/ k. 

Ck is an open set since it is the union of open sets . It is dense on the discrete 
spectrum Sd since it contains Uy for all z/ which are dense in Sd- We have 
Ck+i C Ck- The set £ = flfc -^fc which is a countable intersection of open sets 
Ck dense on 5^, is not an open set but is also dense in Sd (according to Baire 
theorem [25]). When Q ^ C, there exists k such that Q ^ C^, series (49) is 
absolutely convergent and remains bounded at infinite time Thus, the energy 
transferred to the phonon bath remains finite and bounded at all time. This 
set C has zero Lebesgue measure since fi{C) < fi{Ck) < 43/ k for k arbitrarily 
large. Thus, the complementary set 7l\C has full Lebesgue measure. 

By analogy with the theory of Liouville numbers, we now prove that this set C 
contains a (zero measure) subset C (Z C for which series (49) surely diverges so 
that there is (anomalous) energy dissipation. We can also prove that this set C 
does not contain any u^, with probability one if we assume the distribution of 
eigenvalues is uniform enough on Sd- Actually, the probability distribution 
of these eigenvalues is nothing but the measure density of states g{uj)duj. We 
have to assume that this measure is absolutely continuous and moreover that 
the distribution probability of a pair of two eigenvalues uji, and ojyi for v ^ v' 
is obtained as the product of probability distribution of single eigenvalues (i.e. 
the eigenvalues are uncorrected). 

The energy transferred to the phonon bath given by series (49) diverges if its 
sequence of coefficients does not go to zero. This situation surely occurs when 
for any given there is an infinite subsequence uj^^ such that — cjj,,^ | < ka"J^ 
where k is some positive constant. This condition is equivalent to say i7 G £' 
where C! = fliv '^'n ^'n union of open intervals 



^ Note that does not have necessarily to decay exponentially but could also 
decay as a power law with an exponent (3 strictly larger than 3. This is sufficient for 
having 'Y^v''^'^ ^ +oo. Indeed if we assume < < then r^, < ^(g^i^,)/2 so that 
if /3 > 3 this series can be made convergent with an appropriate choice of e > 0. 



43 



,u>N 




(B.3) 



£'jv is a dense open set on Sd and /^tv+i ^'n ^ ■ Thus the countable 

intersection C = f]j^ is a dense set (according again to Baire theorem [25]). 
The Lebesgue measure of this set C is zero because < ^{C'^) < 

4fc Z]i/>Ar('3'!y)^^^ which is the rest of an absolutely convergent series since 
decays faster than l/z/^"*"^, since we assumed an exponential decay. 

We have now to prove that giving an arbitrary uj^,, it could belong to this 
set C only with probability 0. If cu^ G £' we have uj^, G for any N. 
Let us call P the probability for uj^ to belong to C = CI'^'n- l^t us call 
Pn, the probability that belongs to the unions of intervals C'j^. We have 
< P < Pj\f. We choose N > u in order Ui, is not the center of an interval of 
C'^. Then, since we assumed there is no correlation between the probability 
distribution of uJi,/ with that of u,, the probability for uj^ to belong to 
is P/v = Jc'^ g{ijj)duj where g{uj)duj is the density of states. When goes to 
infinity, the Lebesgue measure of C'^ goes to zero. Consequently P < P^ 'iN 
has to be zero. Since the set of is countable, there is no ujy G C with 
probability one. [f]- 



^ This proof should be extendable when the probability correlations between ojy 
and iOyi decay fast enough as |i/ — diverges. 
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C Lemma and its Proof 



Lemma Let us consider a measurable subset £ C TZ of the real line with the 
following property: 

\/u!i e £ and \/u!2 £ £ 

and Vni e Z and 712 & Z (integers positive or negative) 

we have riicui + 712^)2 £ £ (C-l) 

(1) Then, either the Lebesgue measure ii{£) of £ is zero or £ = TZ is the 
whole real line. 

(2) There exists sets £ with property (C.l) with zero Lebesgue measure and 
which are uncountable 

Since the measure of the real hne is infinite, measurable means locally mea- 
surable, i.e. for any finite interval /(a. A) — [a,a-\- A] of width A, the subset 
/(a. A) is measurable. 

There exists many examples of sets having property (C.l): the set of integer 
multiples of a given number, the set of rational numbers. Such sets are count- 
able and thus have zero measure. We shall give also examples of bigger sets 
which arc not the real line, but uncountable. This theorem proves that they 
must have at least zero Lebesgue measure. 

C. 1 Proof of proposition (1) 

Summary of the proof: We assume there exists a measurable set £ with the 
above property which is not the real line and such that there exists an in- 
terval in which it has a non zero Lebesgue measure /^(/(a. A) fl^) 7^ with 
//(/(a, A)n^) < A. Then 

• we prove it has a uniform measure: /x(/(6. A) fl £) only depends on the width 
of the interval A and is proportional to this width (statement 1) 

• for X ^ £ , we consider the sequence of translated sets £n = £ + nx for all 
integers n > and prove that there exists a finite integer p such that £p = 
£q = £. Then, £kp+r = with k and < r < p integers and ^^fl^r' = 
when r ^r' and < r' < p (statement 2) 

• we consider the sequence of translated sets ^„ = £^ + x/p** for n > and 
prove they are all disjoints. Then, ^(/(a. A) nn>i(U-^n) should be infinite 
which is impossible since it is bounded by A. This contradiction proves the 
absence of measurable sets with the required properties and thus proves the 
theorem (statement 3). 
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Let us assume £^ is a measurable set with property (C.l) and with a non 
vanishing measure. 



C.1.1 statement 1 

We prove first that fi{I{a, A) fl^) = ccA where a is some non vanishing con- 
stant between and 1. 

We choose two numbers cui & £ and CU2 & £ such that LO-ijiOx is an irrational 
number. This is always possible because if all 002/ iOi would be rational, then 
the set £/ujq with uoq E £ would be a subset of the rational numbers which is 
countable. Thus, the measure ii{£) would be zero which contradicts our initial 
assumption. 

Thus, we can find Ui E £ and UJ2 £ with UJ2/UJ1 irrational. Then, the set of 
numbers 0Jni,n2 = niUJi + n20J2 G £^ is dense on the real axis (Weyl theorem). 
Thus for an arbitrarily chosen y, we can choose a sequence G £ in that 
dense set {oJni,n2} such that limj^+oo uii = y — x. 

Then, the sequence of translated sets £i = oOi + £ = £ is invariant. Wc have 
obviously by translation {x , A) f] £) = ^{I{x + Ui, A) f]£i) = fi{I{x + 
Ui, A) n^^)- Then, for i — >• +00 n{I{x, A)f]£) = n{I{y, A) fj^)- Consequently, 
n{I{a, A) n £) does not depend on the origin a of the interval /(a. A) but only 
on its width A. 

Then, it comes out by splitting the intervals A = Ai + A2 into arbitrary pieces 
that /x(/(a, Ai + A2)n^) = {a, Ai) []£) + {a, A2) (]£) which readily 
implies that fj.{I{a, A) f]£) — aA is proportional to A independently of the 
origin a of the interval. Nota: It is not trivial that the property f{x + y) — 
f{x) + f{y) for any x and y implies f{x) is linear. However, if f{x) is also 
monotonous (for example increasing) which is obviously true for our function 
IJ>{I{a, Ai + A2) n^) as a function of A, the proof that f{x) = ax is linear is 
straightforward. First, we prove f{x) ~ ax when x is rational and since this 
function is monotone and continuous, we also have f{x) — ax for x irrational. 



C.l. 2 Statement 2 

Let us assume now that £ (Z TZ which implies 3x ^ £. We consider the 
translated set £1 — {£ + x). We have £if]£ = ^. Indeed, if we had a nonempty 
set as intersection , we should find uji E £ and 002 E £ such that uji + x = 
UJ2 therefore, a; G in contradiction with the initial assumption on x. The 
measure of the translated set is identical to that of the original set since it has 
uniform measure: /i(/(a. A) f] £1) = n{I{a — x, A) f]£) = aA 
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Let us now consider the sequence of translated sets £2 = {£ + 2x), £n = {£+ 
nx). If £n f]Sm ^ ^, there exists uji and 002 G S, such that uji+nx = UJ2 + mx. 
Consequently (n — m)x = UJ2 — ^1 & £■ Then and arc identical since 
Em — £'-\- TTix — S + uJi — UJ2 + nx — £ + nx — En- Conversely if (n — m)x ^ S, 
then Enf]Sm^ 0. 

Thus, if for any integer p > 1, px ^ S, the series of sets Sn are all disjoint. The 
union U = Uri=o,oo ^n, is a measurable set and fJ^{I{a, A) ClU) is the sum of the 
measure of the disjoint components J2n f^iH'^j ^) C\^n) = «^ x 00 = +00 
which is divergent. This is impossible, because this measure must remain 
bounded by the measure of the whole interval, A . Consequently, the as- 
sumption px ^ S for any p is wrong. There are integers such that px & S. Let 
us consider the smallest integer po > 1 such that pqx e £. 

Then, Sp^ = £ and if n is written as the integer division n — kpo + r with the 
rest < r < we have £n — Sr and for r 7^ r'we have £rr\Sr' — ^ ■ 



C.1.3 statement 3 

For X E S and integer p > 1, we now consider the infinite sequence = x/p^ 

for n = 0, 1, 2, .... We have Un ^ £ since Hn E £ would implies p'^yn = x E £ 

which contradicts the assumption x ^ £. Otherwise, we have yn — Um ^ £ for 

any n ^ m. It is Vn — Um — x{p^~'^ — l)/p^ ■ Vn — Um ^ £, we would have 

P^iVn — Urn) = x{p"^~^ — 1) E £. Assuming m < n, this is impossible because 
p-m-n _i — i^p_^p_i ^it]^ /;■ = p'n-n-1 ^^^j^ 0<r = p— l<p. Consequently 

the sequence of translated sets J-'n = {£ + y-n) are all disjoint. Therefore, the 
measure /i(J(a, A) nUn>i-^n) is infinite and cannot be smaller than A This 
contradiction proves that £ cannot have nonvanishing measure if it is not the 
whole real axis. 



C. 2 Proof of proposition( 2 ) 

Let us consider a sequence of positive numbers ujn > (n = 0,1,... + 00) 
such that limji^+ooCi;„+i/a;n = 0. This property implies that cUn goes to zero 
with n — > 00 faster than any exponential, that is for any 7 > 0, there exists 
K{^) > such that 

tUn = 1/n! is an example of such a sequence. Obviously, series J2n=o'^n < +00 
is absolutely convergent. 
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Then, let us now consider the set of real numbers £ defined as 

+00 

X & £ : X TJlnOUn (C.2) 

n=0 

where {m„} is any bounded sequence of integers positive or negative that is 
there exists an integer N such that sup„ \mn\ = N < +00. Then series (C.2) 
is absolutely convergent since J2n=o l^nlf^n < NJ2t=o^n < +00. 

Then, we have £ — Uat £n where subsets £jv C £ are defined as 

+00 

X e £n ■ X = ^ rrinUJn 

n=0 

with |m„| < N for all n. We have £n C £n+i 

Thus, assuming x E £ and y E £, there exists A^^^ and Ny such that x G £nx 
and y G £n . Then for any integers and Uy, we have z = n^x + Uyy G 
£\na:\Na:+\ny\Ny which implies z E £. Consequently, £ has property (C.l). 

We now prove that £ has zero Lcbcsguc measure, i.e. = 0. Since £ is a 

countable union of subsets £n, it suffices to prove that ij,{£n) — for any N. 

Giving X — Yl'^=o iTin'-^n £ £n and p an arbitrary positive integer, the partial 
sum Xp — YX^o "'k^'k takes at most (2A^ +1)^ values. The modulus of the 
difference x — Xp is bounded as |a; — = J2n=p l^nl^n ^ NRp where Rp = 
Y^n^pU^n- Then, £i\f is included in the union of (2A^ + 1)^ intervals of width 
2NRp which have a total measure bounded as h{£n) ^ Bp = 2N{2N + lyRp. 

Since for any 7 > 0, there exists -^'(7) such that Un < K{'y)e~^'^, we have 
Rp < K{j)e-^P/{1 - e-^') = K'{^)e-^P. Choosing 7 > ln(2A^ + 1), the upper 
bound Bp < 2NK'{j) {{2N + l)e~'^y goes to zero for p — >• +00 and thus can 
be made infinitely small. Consequently, the Lebesgue measure h{£n) of £n 
must be zero, which implies that £ — Uiv=o^Ar has also Lebesgue measure 
zero. 

For proving that £ is uncountable, it suffices to prove that one of its subset 
is uncountable. This set is defined as x & T : x — J2n=p ^n'^n with m„ = 
for n < p and m„ = ±1 or 0. Since cun+i/ujn has been assumed to go to zero 
for n +00, integer p can be chosen such that cj„+i < cUn/S for any n > p. 
Then for any p < m < n, we have the strict inequalities < a;„ < cu^/S'^"'". 
This choice will ensure that all numbers in T are surely distinct. 

Considering x = I]^^m„c<j„ G and y = I]^^m'^cj„ G where sequences 
{^n} 7^ {^n} different, we define q as the smallest n such that m„ 7^ m'„. 

Then, we have X-y = iinq-mq)Uq + J2n>q{^n-'m'JUn. Since | T,n>q^n^n\ < 

Y.n>q \mn\uJn < (^qJ2n>i 1/3" = (^q/2 and Similarly | T,n>q'm'n^n\ < ^^g/2, we 
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have \x — y — {rUq — m'q)ujq\ = \ Y.n>qij^n — iT^'n)^n\ < ^^q- Siiice rrig and m'^ are 
different, we have |mg — m'^l > 1 which imphes that x and y cannot be equal. 

Since, the set of infinite sequences m„ with m„ = ±1 or for n > p is 
obviously uncountable (it can be mapped to the interval [0, 1] where numbers 
are represented in base 3), JF is uncountable as well as E which proves the last 
statement of proposition (2). 



D Derivation of the integro-differential equation for time-periodic 
solutions 

Searching for periodic solutions with period T — 2'k/^I 

q{t + T;Q) = q{t]Q), u^{t + T;Q) = u^{t]Q) (D.l) 
for all t and all u, allows a representation as Fourier series 

oo oo 

g(t;0)= Yl qnme'''''', u,{t;n)^ ^ 4")(0)e^^-* . (D.2) 

n=— oo n=— oo 

Substituting this into eqs.(16),(17) and ehminating u^J^\Q) leads to 

- {nnrqnm + {V'iqmnm - E "^^l^^^^^l^.^^ - O . (D.3) 
The back Fourier transform of eq.(D.3) yields 

T 

q{t; n) + V'(q(t; Q)) - \^ j G{t - r; Q)q(T; Q)dT = (D.4) 



with 

oo 

G(i;Q) =G(i + T;Q) = ^ G!„(Q)e^""* (D.5) 

n=— oo 

and 
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Note, that ^^"^(Q) is even in n. Accordingly, G{t; Q) is symmetric in t. Prom 
eqs.(D.5) and (D.6) we obtain 



oo 



n— oo 

oo 

= G^^^ + 2j2G^:j,in)cosnnt . (D.7) 



With [35] 



^ cosna; 1 tt cos a(7r — a;) , 

2=^-0 ^ ,0<a:<27r p.8 

— a^^ la^ 1 a sm cktt 



and eq.(D.6) we get from eq.(D.7) 



which is identical to Ga,p{t'i Jl) from eq.(52), taking into account Ga,p{—t] Q) = 
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